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cn . Abstract 

o 

^N I We present the method and compute the strong couphngs of the lowest and first 

orbitahy excited heavy-hght mesons to a soft pion in the static heavy quark limit 
on the lattice. Besides the usual g and g couplings, we were able to make the first 

^ [ computation of the coupling h using the relevant radial distributions. Our results 

are obtained from the simulations of QCD with A^f = 2 light Wilson-Clover quarks, 
combined with the improved static quark action. The hierarchy among couplings 
that emerges from our study is g < g < h. 
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1 Introduction 

A lattice QCD estimate of the pionic couplings to the pair of heavy-light mesons belonging 
to the same or different doublet (s) is important for at least three reasons: (i) They pa- 
rameterize the hadronic matrix element relevant to the kinematically allowed D* — )■ Dvr, 
Dq — )■ Dtt, D'^ — )■ -DqTt, and the similar 5-decay modes when they are kinematically al- 
lowed; (ii) They help us understanding better the viability of various model-approaches 
used in describing the hadronic interactions at low energies, that are still used whenever 
the calculations based on simulations of QCD on the lattice cannot be made. Further- 
more, these couplings are helpful in understanding the dynamics in crossed channels that 
are expected to drive the behavior of the desired D — )■ vr and B ^ n form factors; (iii) 
They enter the expressions derived in heavy meson chiral perturbation theory (HMChPT) 
that are used to guide the chiral extrapolations of the i?-physics quantities computed on 
the lattice. 

The simplest such a quantity is the so called g coupling that essentially represents the 
hadronic matrix element of the light axial current between the meson states belonging 
to the j[ = (l/2)~ doublet. It was computed in the static heavy quark limit on the 
lattice, first in the quenched approximation (A'^^ = 0) [H [2], and later on by including 
Nf = 2 dynamical quark flavors [21 IH El E] • For the propagating charm quark this coupling 
has been computed in refs. [3 E] and it turned out to be considerably larger than in the 
static heavy quark limit. Moreover, the values for ^f-coupling obtained in lattice QCD 
were much larger than the ones obtained by using various QCD sum rules techniques (c.f. 
ref. [9] and references therein) widely used in the -B-physics phenomenology. Finally, the 
result obtained with the propagating charm quark appeared to be consistent with the 
experimentally established vector meson width r(D*) |10j . 

Similar pionic couplings that relate an orbitally excited and the lowest pseudo-scalar 
heavy-light mesons are harder to compute, especially in the static heavy quark limit. To 
our knowledge the only attempt in that direction has been made in ref. [11], where the 
authors computed the B^ — > Bn decay amplitude in the static heavy quark limit directly, 
without relying on the soft pion theorem. 

We remind the reader that the heavy-light meson states belonging to the j[ = (l/2)~ 
doublet are denoted as {B,B*) having J^ = (0~,1~), while the first orbital excitations 
j[ = (1/2)+, are denoted as {Bq,B[), and have quantum numbers J^ = (0"'",1"'"). A 
similar notation is used for the corresponding D-meson states. 

In our recent paper [12] we presented the results of our study of the radial distribu- 
tions of the light quark axial current density sandwiched by the static meson states both 
belonging to either j[ = (l/2)~ or (1/2)+ doublets. That study was a continuation of the 
research made in ref. [13] in which we implemented the improved form of the Wilson line 
on the lattice (static heavy quark propagator) by using the so-called hyper-cubic blocking 
(HYP), a procedure described in detail in refs. [HI dSj. We showed that the accuracy of 
the lattice data can be further improved by applying the HYP blocking on the heavy spec- 
tator quark twice. Later on the radial distributions reported in ref. [12] were compared 
with those obtained in two classes of quark models and a striking qualitative and even 
quantitative agreement has been found [16] . 



The integration of the radial distributions of the hght axial current densities leads to 
the pionic couplings g and g, for the case of heavy-light mesons belonging to j[ = (l/2)~ 
and (1/2)"'", respectively. The fact that we were able to study the radial distributions to 
a good accuracy opened a possibility for computing the axial coupling h, that relates the 
heavy light mesons belonging to two different doublets, such as the one that parameterizes 
the hadronic amplitude in Bq -4 Bit decay, i.e. the scalar to pseudoscalar decay with the 
emission of an S'-wave pion. 

In the following we will describe how the coupling h can be computed on the lattice in 
the static heavy quark limit. We will remind the reader of the difficulties encountered in 
the study of h, and present our solution of the problem. We will then present the numerical 
results from our lattice QCD calculations with Nf = 2 dynamical quark flavors of Wilson 
type. In addition, we will provide the other couplings, g and g, that are then combined 
with /i-coupling in the combined chiral extrapolation. 

2 Methodology 

2.1 Definitions and relevant hadronic matrix elements 

We begin by explaining the method to get the /z-coupling in the soft pion limit. The 
starting point is the hadronic matrix element of the light quark axial current between the 
scalar (-Bg), and the pseudoscalar (B) mesons, 

{B{p')\A,{0)\B*ip)) = ip + p'),A4q') + q,A4q'), (1) 

where A±{q'^) are two Lorentz invariant form factors, functions of the momentum transfer 
q^ = {p — p'Yi and the charged light axial current is A^ = u'^^'^^d, that in the following will 
be denoted as A^ = cpf^'-i^q, reflecting the fact that we take the quarks to be degenerate 
in mass, tti^ = m^ = rUq. After taking the divergence of the axial current in the soft pion 
limit we have, 

hm {B{p'WA,mBl{p)) = {ml, - m^)A+(0) . (2) 

On the other hand, by means of the reduction formula, 
fnm'^ 
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n^{q)B{p')\B*,{p)) = {B{p')\q^A,mB*oip)) , (3) 



so that by deflning 



{7i^{q)Bip')\B*{p))=gB^,B., (4) 

in the same limit, g^ — )■ 0, we have 
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gB*B. = ^^. A+(o). (5) 



The definition which relates the couphng Qb^btt with the couphng h, that appear in the 
chiral Lagrangian [T7], reads 

QBi^Bn = ^rriBmB* ^^ j- , (6) 

and therefore in the static heavy quark mass hmit we can identify, 

h = A+(0) . (7) 

This couphng describes the emission of the S'-wave soft pion off the B^ meson, i.e. 

r(5*^5vr±) = ^|g;|, (8) 



[m|, - {niB - m^)2][m|, - (m^ + m^y] 

Wn\ = ^ ^ , (9) 

and similarly for the charmed heavy-light mesons. 

The main goal of the present paper is therefore the extraction of the form factor Aj^{0). 
To ensure the standard heavy quark effective theory (HQET) normalization of states, 
{Ba{v)\Bi,{v')) = SabS{v — v'), and by taking both in and out states to be at rest, we rewrite 

where we introduced A = ttlb* — i^b, the mass difference between the lowest orbitally 
excited heavy-light meson and its L = counterpart. In the static heavy quark limit 
{mh — 7- oo) we finally have, 

(i?|Ao(0)|5*)^Q^T = A+(A2). (11) 

Notice that the argument of the form factor is non-zero even in the static heavy quark limit 
(A 7^ 0) because the two mesons belong to different doublets of heavy-light meson states. 
In order to go from yl+(A^) to the coupling h = A^{0), an extra step is needed. Before 
we discuss that issue, we will spend some time explaining how the yl+(A^) is computed on 
the lattice. 

2.2 Extraction of A+{A'^) from the correlation functions com- 
puted on the lattice 

On the lattice we first compute the following three-point correlation functions 
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where (. . . )y denotes the average over independent gauge field configurations, the inter- 
polating operators for the pseudoscalar and scalar heavy-light mesons are respectively 
P = h'-f^q and S = hlq, with h{x) and q{x) being the static heavy and the light quark fields 
respectively. The purpose of computing the elastic correlation functions, Cpvp,svs{ty,tx) 
with Vq = qjoq, is to eliminate the interpolating source operators. Those elastic matrix 
elements are fixed by the electric charge conservation and, after accounting for the ap- 
propriate renormalization constant, the relevant form factors are simply equal to one. To 
extract A+(A^) we adopt the strategy of double-ratios and compute 
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{B\Vo\B) {B*\Vo\B*) 

where an extra "— " sign accounts for the fact that the matrix elements for the pion emission 
and the pion absorption have opposite signs. In terms of quark propagators 
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where we used the 75-hermiticity of the quark propagator in the euclidean space, SUx, y) = 
'l5<Sqiy,x)j5. In addition, we wrote the static quark propagator as (1 -|- 7o)/2 x W^, with 
the Wilson line being 
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The latter is merely obtained from the discretized static heavy quark action [18 
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where for the time component of the link variable, Uq , we use several improved schemes 
discussed in our previous paper |T2]. Improvement of the discretized Wilson line is essential 
as it ensures the exponential improvement of the signal to noise ratio in the correlation 
functions with respect to what is obtained by using the simple product of link variables 
Similarly, other correlation functions are computed as 
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In the following we drop the dependence on ty as this time is kept fixed to one or several 
values that are sufficiently large so that the lowest lying state would be isolated. In 
practice we actually check that by reversing the roles of tx and ty our results remain 
stable. The use of improved static heavy quark actions is essential because we also need 
the signal to be good enough at large time separations between the source operators. 
However some efficient smearing of the sources is welcome in order to have both sources 
sufficiently far away from the light quark current inserted between them, and therefore 
that the desired matrix element can be isolated. We use the smearing proposed in ref. [20] , 
that we already discussed in our previous papers [1], which consists in replacing the source 
operators h{x)rq{x) — > h{x)rq^{x), where F is a Dirac matrix, and 
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1/2)^. After inspection we find the smearing to be efficient for 



i? = 3 and -Rmax = 5. To make sure that the values we obtain are correct, we also tried the 
standard method of extracting the form factors from the three point functions, by dividing 
out the source operators obtained from the fit to two-point correlation functions computed 
on the lattice, i.e., 
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We made that exercise by including either one or two exponentials on the right hand side 
in order to check the stability of our results. With all these checks we were able to find 
the stability window in time, that indeed leaves the form factor ^^(A^) unchanged. 

2.3 Going to q^ ^ 

As we already mentioned, the information we obtain from the study of the three point 
correlation functions with both source operators at rest is the form factor yl+(A^), while 
we actually need A+{0). Going to g^ = would require a momentum injection to one of the 
source operators, which is particularly complicated to do in the case of static heavy-light 
mesons. Furthermore, tuning the three-momentum to \p — ff\ = A is very difficult too. 
This is where our experience with radial distributions of various light quark current matrix 
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Figure 1: Shapes of the distribution fpAs{f), as extracted horn the three and two point correlation 
functions computed on the lattice as indicated in eq. Il21\) . lUustration is provided with one value of the 
sea quark mass corresponding to Kgea = 0.1374 from Set-1 lattices listed in tab. [TJ and with the HYP^ 
static heavy quark propagator. All data points are given in lattice units. 



elements becomes useful. We follow the same procedure employed in ref. [12] and compute 
the radial distribution of the correlation function CpAs^ namely 



CpAsity, t,; r) = (Y, Hxh5q{x) g(0 + f)7o75g(0 + F) (/i(y)Ig(y))t)^ 
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that leads to 
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for the sufficiently large time separations among operators. Note that we place one source 
operator at t^; < and the other at ty > 0, while the light quark axial current, Aq = g7o75(?, 
is fixed at the origin of the lattice. The shape of the function fpAsi^) for one representative 
value of the light quark mass is shown in fig. [H U 
We checked that the integral, 
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indeed reproduces the form factor value obtained from the double ratio (TT^ . even if with 
slightly larger error bars. To get the form factor at q^ ^ A^ we then give the momentum 



"'^Note that in our study we always consider the sea and the valence light quarks to be mass degenerate. 



q = (0, 0, g^), by taking the Fourier transform, 

//•OO /"TT 

df fpAsir) e"^ = 271 dr r^fpAsir) e''^^^""''^ sin 9d9 
Jo Jo 

= Att rdrr'^^^^^fpAsir). (23) 

Jo QzT 

Obviously, by taking qz = A = S — £, that we easily extract from the study of the time 
dependence of two point correlation functions, as indicated in eq. (TTOl) . we get the desired 
form factor A^{0), namely 

^+(0) = 4vr r dr r''^^fpAs{r) . (24) 

Jo Ar 

Since the mass dependence in renormalization constants with improved Wilson action on 
the lattice cancel to a large extent in the double ratio ( fT3|) . we prefer to determine the form 
factor A+(A) in that way. Indeed, from the above eqs. (I23II24I) . we can get the correction, 

_ A+{0) _ Jo Ar 

' dr r fpAs{r) 



and therefore our value for the desired pionic coupling h will be obtained via, 

hg = A+{Al)xR^^, (26) 

where in the last line we added an index q to distinguish hg from the true coupling h that 
is obtained in the chiral limit. 

2.4 g and g couplings on the lattice 

For a consistent chiral extrapolation of hg, we also need the couplings g and g, that we 
discussed in our previous papers [H |12]. With respect to ref. [12], here we take the light 
quark current to be at the origin of the lattice, which helps avoiding the fact that the ratio of 
correlation functions becomes a sum of all possible diagonal matrix elements. Furthermore 
we extract the couplings from the ratios of three point correlation functions instead of three 
to two point functions. In short, the couplings g and g are defined via 

{BiO)\q^a5q\B*{e,6))=e,g, 

{B*,{0)\q-f,^,q\B[{e',6))=e',g, (27) 

where e^ and e' are the polarization vectors of the vector and axial- vector heavy-light me- 
son, respectively. To extract the above couplings from the ratios of three-point correlation 
functions we also computed: 
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which when combined with those hsted in eq. ( 1T2|) lead to 

R (f f \ — — ^PAv{ty:tx) {B\Ai\B*) _ ^ 

Zy Cpvp{ty,tx) {B\Vq\B) 
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J^gy^yyi^x) '7 ^ (+ +\ /R*lT/lR*\ ■^'^ ' ^ ' 

where we again added an index q to distinguish cjq and 'cjq from the true pionic couphng that 
are defined in the soft pion hmit. To our knowledge the computation of these couplings 
from the ratios of three-point correlation functions alone has not been attempted elsewhere. 

3 Lattice computation of the pionic couplings 

In this section we present our numerical results. We use three sets of publicly available 
gauge field configurations, all obtained by using the Wilson actions. In particular, the 
first set of data is obtained by using the Iwasaki gauge action, while for the other two 
sets the standard Wilson plaquette action has been used. As for the discretized Dirac 
part of the QCD action, in all three sets of configurations the Wilson-Clover action was 
implemented to ensure the elimination of 0{a) discretization effects from the lattice action. 
More specifically, in the first set the Clover coefficient has been set to its value obtained 
by using the one-loop boosted perturbation theory {csw = 1-47), while in the other two 
ensembles the non-perturbatively determined Clover coefficient was used in simulations 
{csw = 1.919, and 1.823, for (3 = 5.29 and 5.40, respectively). Obviously, in computation 
of the light quark propagators we used the same values of csw used to generate the gauge 
field configurations with A^f = 2 fiavors of mass degenerate light quarks. From the set 
of publicly available gauge field configurations we chose those that are separated by 20 
unit-length HMC trajectories. 

More details concerning the simulations from which the gauge field configurations used 
in this work have been obtained can be found in refs. [21], [22]. In tab. [1] we summarize the 
main details concerning the gauge ensembles used in this work. To situate the values of the 
lattice spacings, one can use Tq = 0.47 fm [23], which then translates to a"^ ^ 2.0, 2.6 and 
3.1 GeV, for our three sets of gauge field configurations respectively. It should be stressed 
that the value of tq is unphysical and is still a subject of uncertainty. For example, from 
simulations with twisted mass QCD on the lattice and by using the physical pion decay 
constant one obtains tq = 0.44(1) fm [21]. For our final results, the uncertainty in tq is 
unimportant as our main goal is to obtain the dimensionless pionic couplings in the chiral 
limit. 

In tab. [T]we also give the values oi Z = ZA(gQ)/Zv{gQ) that is used to extract A+(Ag), g, 
and g, from Rh, Rg, and Rg respectively. Besides the cancellation of the exponential terms 
these ratios are convenient because the sizable mass corrections to the renormalization 
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m^ [GeV] 


ref. 


Z (NP) 


Z (TI) 


Iwasaki/ Clover 


2.1 [4.70(2); 


0.1357 


1.025(5) 


|21| 


1.03(2) 


1.017 


Set 1 




0.1367 
0.1374 
0.1382 


0.888(4) 
0.761(4) 
0.563(4) 








WP/Clover 


5.29 [6.20(3); 


0.1355 


0.858(6) 
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1.038 


1.029 


Set 2 




0.1359 
0.1362 


0.643(7) 
0.406(6) 








WP/Clover 


5.40 [6.95(4) 


0.1356 


0.969(6) 


II22I 


1.033 


1.028 


Set 3 




0.1361 
0.13625 


0.672(7) 
0.572(8) 









Table 1: Gauge Held configurations with Nf = 2 dynamical Wilson-Clover quarks, where the hopping 
parameters Ksea is specified together with the corresponding pion mass. "WP" stands for the Wilson- 
Plaquette gauge action, and Clover stands for the O{a)-improved Wilson quark action. All lattice volumes 
are 24'^ x 48. More information on each set of configurations can be found in the quoted references. The 
ratio of renormalization constants is denoted by Z ~ Z Aioo) / Zy {go) . "NP" stands for non-perturbatively 
determined value, that is compared to a tadpole improved perturbative one ("TI"). 



constants cancel to a large extent in the ratio. More specifically, for the consistent 0{a)- 
improvement of the Wilson quark action, the renormalization constants are improved as: 
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where Zy ji^^g^) stands for the value obtained in the chiral limit, and 6y,A(fi'o) is the counter- 
term coefficient chosen to eliminate the the 0{amq) effects when working with the vector 
and axial- vector current respectively. Knowing that by ~ b^, the effect of 0{amq) cor- 
rections to Z = ZA{gQ)/Zv{go) becomes negligible. Furthermore the non-perturbatively 
determined value of Z in the chiral limit has been found to be very close to its value 
predicted by using the tadpole improved one-loop perturbation theory [23]. This can be 
appreciated from the numbers in tab. [T] where the non-perturbatively computed Z [2S] are 
less than 1% larger than the corresponding tadpole improved perturbative estimate [251126] . 
In the first set instead the non-perturbative value for the ratio was deduced from the Ward 
identities [1], which turns out to be very close to the perturbative result [21]. In what 
follows, for the first set we will use the perturbative result (TI), while for the other two 
the non-perturbative values will be used. 

Our main results will be obtained from Set 1 of the lattice data for which we were able 
to compute all three couplings. With the other two ensembles the computation of gg was 
unstable, but the computation of cjg and hg is good which is why we prefer to separate the 
discussion of our results. 



3.1 Discussion of the numerical results 

We summarize our numerical results in tab. [2] that we comment in the following. The mass 
difference between the lowest lying (1/2)+ and (l/2)~ states, A^ = £^g — £g, is obtained 
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0.1357 


0.269(4) 


0.639(5) 


-0.040(10) 


0.809(24) 


0.889(26) 


0.719(29) 


0.1367 


0.273(4) 


0.621(3) 
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0.832(14) 


0.896(15) 
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-0.103(12) 


0.889(26) 


0.842(24) 


0.751(30) 


0.1382 


0.272(6) 


0.529(19) 


-0.171(19) 


0.968(39) 


0.808(32) 


0.783(44) 


0.1355 


0.232(8) 


0.552(19) 


-0.094(10) 


0.981(38) 


0.829(32) 


0.814(45) 


0.1359 


0.210(5) 


0.514(17) 


-0.143(17) 


0.892(31) 


0.797(28) 


0.711(35) 


0.1362 


0.223(12) 


0.454(23) 


-0.171(20) 


0.919(56) 


0.772(47) 


0.710(61) 


0.1356 


0.216(6) 


0.573(15) 


-0.110(11) 


0.978(40) 


0.840(35) 


0.822(48) 


0.1361 


0.233(7) 


0.519(24) 


-0.155(11) 


0.903(32) 


0.822(29) 


0.743(37) 


0.13625 


0.207(8) 


0.483(19) 


-0.164(8) 


0.868(23) 


0.810(21) 


0.703(26) 



Table 2: Direct numerical results extracted from the correlation functions calculated on all of the en- 
sembles of the lattices with parameters listed in tab. [H Results are presented in the same order as in 
tab.H 

from the fit of the form given in eq. ( TT9|) to the two-point correlation functions computed 
on the lattice. All results presented in this section are obtained by using the HYP-2^ static 
quark action discussed in our previous paper [12]. With other forms of the static quark 
actions we obtain fully compatible results with those listed in tab. [2J The benefit of using 
the improved static quark actions is that the signal remains good at larger time separations 
so that we could check that with either local or smeared sources we could obtain results 
that are compatible. Importantly the smearing helps to isolate the lowest lying state at 
shorter time separation between the source operators which is important for the studies of 
the three point functions. We make several tests by fitting each of the three point functions 
to its hadronic decomposition, 
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(31) 
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with either i,j = l or i,j = 2, where the coupling to smeared source operators (2f)j and 

{Z^)j are obtained from the fits to two point correlation functions. In this way we could 
determine the stability interval where the fit to our data lead to a stable value of /1+(A^). 
That is then compared with what we obtain by using the double ratio. Another good 
test for the goodness of the fitting interval is to compare the effective mass plots. We fix 
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—tx = ty = t and since the light quark axial current is fixed at the origin of the lattice, the 
effective mass from three point functions CpAsit^ ~t) = CpAsi^)'- 

'CpAs(t + 1) 



should be compatible with the sum of what we obtain from the two point functions 
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Since both interpolating field operators are pulled away from each other in eq. ( I32l) by one 
time unit, that makes twice larger distances of the source operators in time, than it is the 
case with two point correlation functions obtained from eq. ( !33|) where one source operator 
is kept fixed at t^ = 0. After accounting for that fact we checked that for each of our 
lattices and for every three point correlation functions, the effective mass of three point 
functions is consistent with what is obtained from the two point functions. An illustration 
is provided in fig. [2l which is a typical situation with the lattice data. As expected the 
signal for the excited [(1/2)+] state is less good than the one for the ground static heavy- 
light state [(1/2)"]. The values of A+{Aq), c/g, and 'gg are then obtained from the ratios flT3|) 
and ( l29l) . As we mentioned above, besides the double ratios we also used the ratios with 
two point correlation functions similar to the one in eq. ( 12T|) to check for the consistency 
of the obtained results. This also helped us determine the sign of the form factor A+{Aq). 
After a detailed inspection, we fit our data between —t^ = ty = t & [4,7] for A+{Aq), 
between t G [3, 9] for cjq, and t G [3, 5] for Qg. Finally, we should emphasize once again that 
all the results presented here are fully unquenched, i.e. the valence quark mass is equal to 
the sea quark mass. 

Concerning the distribution fpAs{f), we obtain it by fitting all our data as indicated 
in eq. f l2T|) on the same interval —t^ = ty = t & [4, 7], and for each value of r. The result is 
illustrated in fig. [1] where we show both fpAs{f)i and r'^fpAsi^)- ^^ practice we also check 
in each situation that the sum over all lattice points indeed reproduces the result obtained 
by using the radial distribution. We then checked that 

/■CXD 

A+iA')=A7T drr^fpAsir), (34) 

Jo 

coincides with the value obtained from the double ratio ( !T3|l . even if with a slightly larger 
error. □ Like we discussed in ref. [12], one might worry whether or not Air obtained from 
the angular integration might be spoiled by the fact that we are working on the cubic 
lattice which, of course, is not spherically symmetric. Thanks to the fact that the radial 
distribution fpAs{f) is a fast decreasing function and becomes compatible with zero before 
reaching the lattice boundaries, the effect of the cubic lattice geometry is irrelevant. This 
is the case for the range of pion masses considered in this work. To illustrate that point 
we show that 

hk = 47r V r'^^^'^fpAsir) ^ 47r / dr r'^^^'^ f p As{r) , (35) 



r=0 



^In doing so we of course accounted for the 0{amq) effect in the renormahzation constant. 
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Figure 2: Effective mass plots: £q '{t) and Eq ' {t) refer to the (l/2)~ and (1/2)+ states respectively, 
both obtained from the two-point correlation functions computed on the lattice as indicated in eq. i33\) : 
{£q + £q){t) is obtained from the three-point correlation functions, c.f. eq. ij32)) . Dashed lines correspond 



to the central values of the static heavy-light binding energy: Eq , Eq , and Eq + Eq . Illustration is 



(2) ?r(2) 



7(2) 



:^{2) 



provided for the lattices of Set 1, with Kg. 



0.1374, and by using the HYP-2 static heavy quark 



propagator. All data points are given in lattice units. 



13 








0.0 <^ 




Figure 3: Summing over all distances: Plots show the saturation of the sum at larger distances of Iq that 
corresponds to —hq/ZA, a,nd of I2 that is related to the second moment of the radial distribution. Dashed 
lines are the final results for Iq,2 computed according to eq. i|35|) . Illustration is provided for the lattices 
of Set 1, with Ksoa = f 



0.1374, and by using the HYP-2^ static heavy quark propagator. Distances are 



given in lattice units. 



for k = 0, 1, i.e. that for larger distances the sum becomes saturated and is a good 
approximation to the continuum integral. I2 above is related to the second moment of the 
radial distribution fpAs{f) that we will discuss later on. Illustration for the saturation of 
both Iq and I2 on the specific example of the lattice data is shown in fig. [3l 

Finally, we should comment on the values of hq given in tab. [2] that is obtained in a 
way discussed in subsec. 12.31 namely as: 



hq = A^iAl) R. = A^iAl) 



^,,2sin(A,r) 



Aqr 



fpAs{r) 



dr r'^fpAsir) 



(36) 



This is quite a unique opportunity in hadronic physics that we can check for the moment 
expansion. To that end we expand the integrand in eq. ( 1391) for small QzT and then by 
choosing q^ 



Ag we have 



A+(0) = 47r 



/ dr r'^fpAsij) - -Aj / dr r^fpAs{r) 



120 



A+(A 



A^ A^ 

6 ^ / 120 ^ / 



1 r°° 

-Ajy^ drr'fpAsir) + ... 

(37) 



and after keeping the first couple of moments we see that 



A2 A^ 

R^ ^R^ = l-^/r^) + ^L/r^) 



(3^ 
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where the moments are 



/■oo 



/■oo 

/ t/r r'^fpAsir) 
Jo 



(39) 



We computed the moments (r^) and (r"^) and observe that the first moment alone saturates 
i?A, to about 80%, while R^ roughly coincides with -Ra,- The table of results for the 
moments, as well as a comparison R^/R^ is provided in Appendix of this paper. 



3.2 Chiral extrapolation 

To get the physically interesting couplings h, g and g we need to perform a chiral extrapo- 
lation of the results listed in tab. [21 Since our results are obtained for the pion masses not 
too close to the chiral limit and exhibit a rather clean linear dependence in rriq oc m^, we 
first attempt a linear chiral extrapolation. In particular. 



9q 



g{l + Cgml) 



9q 



g{l + CgTn 



I) 



h[l + Chjnl) 



(40) 



Another possibility is to use the expressions derived in HMChPT [27j, in which the effects 
of the nearest orbital excitations has been taken into account. Those corrections read: 



9q = 9 



9q = 9 



hq = h 



w 



(47r/)2 
A-g" 



2 1 2 

m^ log m^ 

7T <^ 7T 



h' 



mz 



(4vr/)2 8A 



ml log ml + 



h' 



m„ 



(4vr/)2 8A2 



3 + - ) m^ log m^ + Cam,'^ 



-^3"^ 



2 1 2 I ~ 2 

m^ log m^ + Cgm^ 



1- 



^^g' + ^g'-2gg 



4 



(4vr/)^ 



2 1 2 

m^ log m^ — 



?n; 



(4vr/)2 2A2 



^m^logm^ + c,,m^ 



(41) 



where, again, on the left hand side are the quantities we computed on the lattice (each 
carrying an index "g") and on the right hand side are the desired quantities /i, g and 7j. 
The terms oc m^/A^ come from the inclusion of the heavy-light states of opposite parity 
in the chiral loops. Omitting those terms is equivalent to assuming that m-j^ <^ A, which 
is the approximation needed to recover the usual HMChPT formulas for g, and 'g. To 
appreciate the difference, we will fit our data to the form in which the terms oc m^/A^ are 
omitted and the corresponding results will be labelled as "chi-1". The results obtained by 
using the full formulas f HT]) will be denoted as "chi-2" . The fits are illustrated in fig. IH 
and the results are given in tab. [31 

We did not combine all the data to make the final extrapolation. Instead we prefer 
to give our results separately for the data in which Iwasaki gauge action has been used 
(Set 1), and those that are obtained from the field configurations generated by using the 
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Figure 4; Chiral extrapolations: Dashed lines denote the linear extrapolations, and the full curves 
correspond to the results obtained by fitting the data to the expressions derived in HMChPT, c.f. eq. ((47)) . 
In the lower plot we combined the results from Set 2 (lighter) and Set 3 (thicker symbols). Results of the 
fit are given in tab. [3 
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Wilson plaquette gauge action (Set 2 and 3). Although the results are compatible, within 
the error bars, we note that hg values in the Set 1 remains quite flat when the light quark 
mass is varied, while those obtained in Set 2 and 3 exhibit a positive slope in the light 
quark mass. As our final results we therefore present h, g, g separately for Set 1, and for 



o 
CO 



Extrapolation 


h g 


9 


Ch 


Cg 


Cg 


lin. [eq. ()40|)] 


0.80(2) 0.53(3) 


-0.20(3) 


-0.10(2) 


0.22(8) 


-0.80(6) 


chi-1 [eq. (l41l)] 


0.69(2) 0.43(1) 


-0.19(2) 


0.05(4) 


0.51(5) 


-0.80(7) 


chi-2 [eq. (l41j)] 


0.63(3) 0.41(1) 


-0.19(2) 


0.13(8) 


0.59(3) 


-0.80(6) 



CO 

a 

CM 

o 

CO 



Extrapolation 


h g 


9 


Ch Cg 


Cg 


lin. [eq. (|40j)] 


0.64(2) 0.44(2) 


-0.20(2) 


0.31(6) 0.39(9) 


-0.58(10) 


chi-1 [eq. (|41j)] 


0.58(2) 0.40(1) 


-0.22(2) 


0.47(9) 0.52(6) 


-0.57(11) 


chi-2 [eq. (|41j)] 


0.56(2) 0.40(1) 


-0.19(2) 


0.47(10) 0.53(5) 


-0.57(10) 



Table 3: Results of the extrapolations of our results to the chiral limit. 



Set 2&3, namely: 
Set 1 



h = 0.69(2) K9 , g = 0.43(1) K) , g = -0.19(2) {tl) 



Set 2 & 3 : h = 0.58(2) (+6) , g = 0.40(1) (+^) , g = -0.22(2) (tl) 



(42) 



where as a central we take the values obtained from the extrapolation "chi-1", and the 
difference with the results obtained through the linear extrapolation and by using the full 
formulas from eq. (HTl) are used as the systematic error estimate. We believe it is better to 
leave these errors as such, and let the lattice data decide the preferred extrapolation curve 
by investing into simulations at ever lower light quark masses. We see for example from 
the data in Set 2 and 3 that the results obtained from the chiral extrapolations "chi-1" 
and "chi-2" are practically indiscernible, which is a consequence of the fact that the pion 
masses explored through these simulations were lower than those in the Set 1, where the 
difference remains pronounced for the coupling h. 

The emerging pattern from our results is quite clear, g < g < h. Knowing that 
h is rather large, its impact on the chiral extrapolation of phenomenologically relevant 
quantities computed on the lattice should be examined carefully. We will return to that 
question in a separate paper [28]. Here we will simply comment on the decay width for 
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the ^-wave pion emission in the decay of the scalar Dq and i?Q-niesons . If we neglect the 
l/rricfi corrections to /i, then by using eq. (|8]) we obtain 

V{Dl -^ Dn^) = [0.26(13)set i,0.18(6)set 2&3] GeV , 

T{B* -^ Bn^) = [0.22(7)set i, 0.13(6)set 2&3] GeV , (43) 

where we took the physical masses available from PDG [32], except for the unmeasured 
scalar Bq whose mass we get by imposing rriB* — ttlb = ttlbi — ms*-, to get ms* = 
5.68 GeV. A larger error in the charm sector comes from the larger error on the mass 
niD* = 2.32(3) GeV, and the results can be compared to the PDG average of the experi- 
mentally established width of the scalar meson, r(DQ°) = 267(40) MeV [32]. In the case 
of i?Q-decay our results agree with ref. [TTj . 

Before concluding we should also comment on the comparison between our results (H2|) 
and the values obtained by using the light cone QCD sum rules (LCSR) [29] summarized 
in ref. [30] in which also h has been computed: U 

LCSR : /i = 0.60 ±0.13, ^ = 0.17 ±0.06, 5^ = -0.10 ± 0.02 . (44) 



Other QCD sum rule calculations lead to similar results [31] from which we see that the 
LCSR values for g and g are by about a factor of two smaller than those obtained on 
the lattice. This problem persists in the case of the charmed heavy quark, as discussed 
in ref. [9], while the lattice QCD results of go'Dn are consistent with the one extracted 
from the experimentally measured r{D*^) [10]. In contrast to that situation, the value of 
the coupling h obtained from LCSR agrees quite well with the lattice results. Finally, the 
results obtained in two classes of quark models discussed in ref. [16] also agree with the 
couplings presented here. 

4 Summary and conclusions 

In this paper we discussed the computation of the coupling of two heavy-light mesons to a 
soft pion in the static heavy quark limit. We focussed on the heavy- light mesons belonging 
to the j^ = (1/2)- [L = 0], and j^ = (1/2)+ [L = 1] doublets. To compute the coupling 
that relates the heavy-light mesons belonging to two different doublets, h, one encounters 
a difficulty in getting to the point in which the four-momentum of the associated pion is 
g^ — 7- 0. We solve that difficulty by computing the radial distribution of the matrix element 
of the light quark axial current between the scalar and pseudoscalar mesons, which then 
allowed us to compute the suppression factor needed to convert the form factor yl+(Ag) 
to hq = A+(0). In addition to this coupling we also computed the axial couphngs that 
involve the heavy-light mesons belonging to the same doublet, namely g for (1/2)" states, 
and g for the (1/2)+ ones. With all three couplings we were able to check on the chiral 
extrapolation in which the impact of the nearest excitations has been taken into account. 
All results presented in this paper are obtained from the correlation functions computed 



■^ Since the overall sign of the couplings discussed here is convention dependent, we changed the signs 
of the results reported in ref. [30] to be consistent with ours. 
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on the publicly available gauge field configurations obtained with A'^f = 2 dynamical light 
quark flavors and by using the Wilson-Clover quark action and by using the Iwasaki (Set 1) 
and Wilson plaquette gauge actions (Set 2 & 3). We implemented the improvement of the 
spectator static quark propagator (Wilson line), and from the data in which the Iwasaki 
action has been used we obtain: 

Set 1 : /i = 0.71(2)(9), ^ = 0.47(1)(6), 5^ = -0.19(2)(1) , (45) 

where we symmetrize the systematic error that comes from the chiral extrapolation. From 
the data obtained by using the Wilson plaquette gauge action, the results we get after 
combining the values of the couplings computed at two fine lattices, read: 

Set2&3 : /i = 0.60(2) (4) , ^ = 0.42(1)(2), 5^ = -0.23(2)(1) . (46) 

Impact of the results of the present paper on the chiral extrapolation of the phenomenolog- 
ically interesting quantities involving heavy-light mesons as computed on the lattice, will 
be addressed in a separate paper [28] . 
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Appendix 

As we mentioned in the text, thanks to the fact that we were able to compute the radial 
distribution of the desired matrix element, we were also able to check on the moment 
expansion. To do se we used eqs. ( I37t [38 |[39|) . and obtained the results listed in tab. HI All 
results are given in lattice units. 



p 


Kq 


A. 


-Ra 


{r'). 


{r'). 


{Ri/Rd.)q^o\ 


2.1 


0.1357 


0.269(4) 


0.889(26) 


17(1) 


606(26) 


93(4) 


Set 1 


0.1367 


0.273(4) 


0.896(15) 


20(1) 


891(27) 


90(3) 




0.1374 


0.262(6) 


0.842(24) 


22(1) 


1127(65) 


93(5) 




0.1382 


0.272(6) 


0.808(32) 


28(1) 


1835(108) 


92(7) 


5.29 


0.1355 


0.232(8) 


0.829(32) 


30(1) 


1976(100) 


98(8) 


Set 2 


0.1359 


0.210(5) 


0.797(28) 


33(2) 


2542(185) 


99(6) 




0.1362 


0.223(12) 


0.772(47) 


33(2) 


2502(177) 


101(13) 


5.40 


0.1356 


0.216(6) 


0.840(35) 


31(1) 


2033(100) 


95(7) 


Set 3 


0.1361 


0.233(7) 


0.822(29) 


35(1) 


2609(149) 


91(8) 




0.13625 


0.207(8) 


0.810(21) 


39(2) 


3249(180) 


95(10) 



Table 4: The values of the first couple of moments as obtained from our lattices are given in lattice 
units. The comparison of the full expression i?A and the one obtained by combining the moments i?4, as 
indicated in eq. psp , is presented in the last column. For convenience we also give the values of the orbital 
splitting, i.e. A^ — £q ~ £q which is practically the mass difference between the scalar and pseudoscalar 
static-light mesons. Three Sets of data corresponds to parameters listed in tab. [TJ 
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